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Abstract 

We study the energy dependence of the transverse momentum (px) spectra for 
charged pions, protons and anti-protons for Au+Au collisions at ^/snn = 62.4 and 
200 GeV. Data are presented at mid-rapidity (|y| < 0.5) for 0.2 < px < 12 GeV/c. 
In the intermediate px region (2 < px < 6 GeV/c), the nuclear modification factor 
is higher at 62.4 GeV than at 200 GeV, while at higher p T (p T > 7 GeV/c) the 
modification is similar for both energies. The p/ir + and p/ir~ ratios for central 
collisions at ^/snn = 62.4 GeV peak at pr — 2 GeV/c. In the px range where 
recombination is expected to dominate, the p/tt + ratios at 62.4 GeV are larger 
than at 200 GeV, while the p/ir~ ratios are smaller. For px > 2 GeV/c, the p/ir~ 
ratios at the two beam energies are independent of px and centrality indicating 
that the dependence of the p/tt~ ratio on px does not change between 62.4 and 200 
GeV. These findings challenge various models incorporating jet quenching and/or 
constituent quark coalescence. 



Key words: Particle production, recombination, fragmentation, jet quenching, 
nuclear modification factor and particle ratios. 



1 Introduction 



Experiments at the Relativistic Heavy Ion Collider (RHIC) [1] at Brookhaven 
National Laboratory have shown that hadron production at high transverse 
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momentum pr (pr > 6 GeV/c) is suppressed for central Au+Au collisions rel- 
ative to nucleon-nucleon collisions or peripheral Au+Au collisions [2,3]. This 
suppression is thought to be related to jet quenching in dense partonic mat- 
ter [4]. At intermediate p T (2 < p T < 6 GeV/c), in central collisions, the 
baryon to meson ratio is higher than in peripheral collisions [5,6]. This feature 
may be due to hadronization through the recombination of quarks [7]. 

The energy loss by energetic partons traversing the dense medium formed in 
high-energy heavy-ion collisions is predicted to be proportional to both the 
initial gluon density [8] and the lifetime of the dense matter [9]. The energy de- 
pendence of the nuclear modification factor (NMF, defined later) significantly 
constrains parameters in theoretical model calculations. The quantitatively 
large suppression pattern observed at high p T , for both light hadrons and 
those involving heavy quarks [10], has renewed interest in the mechanism of 
energy loss, namely, the relative contribution of radiative and collisional forms. 
The dominance of one over the other depends upon pp and energy [11,12]. Re- 
cently, for a given beam energy a universal dependence of high pp NMF on the 
number of participating nucleons (iV" part ) was proposed as a signature of radia- 
tive mechanisms being the dominant energy loss processes [13]. On the other 
hand, it was suggested that radiative energy loss will break xt (= 2 Pt/V^nn) 
scaling [14]. Thus, a detailed study of the energy, pp, and N pSuTt dependence 
of identified hadron production and hadron scaling properties is needed to 
continue the investigation of energy loss mechanisms. 

In this letter we report the results of such a study performed using iden- 
tified charged pions, protons, and anti-protons for rapidities \y\ < 0.5 and 
Pr < 12 GeV/c for Au+Au at v /si^ = 62.4 and 200 GeV. The data were 
taken by the STAR experiment at RHIC [15]. 

Identified particle pr spectra at different beam energies will also enable the 
study of the effects of the energy dependence of parton energy loss and initial 
jet production on the produced hadron pr spectra. At high px (pr ^ 6 GeV/ c), 
pions are expected to originate dominantly from quark jets at ^/snn = 62.4 
GeV, while both gluon and quark jets contribute substantially to pion pro- 
duction in the same pr region at ^/snn = 200 GeV [16,17]. Therefore, a factor 
of ~ 3 difference in xt (for a given pr) at the two beam energies may allow 
the study of the difference in energy loss mechanisms for quarks and gluons. 
This difference in energy loss is due to the non-Abelian feature of color charge 
dependence of parton energy loss [18,19]. Alternatively, as p production is 
dominantly from gluon jets, the pip) /it ratios are sensitive to quark and gluon 
jet production in heavy-ion collisions [20,21]. Baryon production relative to 
meson production is also sensitive to baryon transport and energy densities. 
The energy dependence of the baryon-to-meson ratio will address the specific 
prediction of the quark coalescence models of a higher baryon-to-meson ratio 
at ^snn = 62.4 GeV compared to 200 GeV in the intermediate pr range [22]. 
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2 Experiment and Analysis 



The data presented here were taken at RHIC in 2004 using STAR's [15] Time 
Projection Chamber (TPC) [23] and a prototype Time-Of-Flight (TOF) [24] 
detector. The TPC magnetic field was 0.5 Tesla. Using a minimally biased 
trigger (MB), 7.4xl0 6 and 1.4xl0 7 Au+Au events at v /s^ = 62.4 and 200 
GeV, respectively, were analyzed. 1.5xl0 7 200 GeV Au+Au events from a 
central trigger were also analyzed, which corresponds to the top 12% of the 
total cross section [21]. The identified particle spectra for Au+Au collisions at 
200 GeV are presented in Ref. [21]. Centrality selection at 62.4 GeV utilized 
the uncorrected charged particle multiplicity for pseudorapidities | r] \ < 0.5, 
measured by the TPC [21,25]. Ionization energy loss of charged particles in the 
TPC was used to identify n ± , p and p within | rj \ < 0.5 and full azimuth, for 
p T < 1.1 GeV/c and 2.5 < p T < 12 GeV/c. Detailed descriptions of TPC par- 
ticle identification techniques for the low px range (0.2 < px < 2.5 GeV/c) 
can be found in Ref. [26]. For px > 2.5 GeV/c, the relativistic rise of ioniza- 
tion energy loss was used to identify the 7r ± , p and p [21,27]. The TOF data 
allowed pion and proton identification up to px ~ 3 GeV/c for —l<i]<0 and 
A$ < tt/30 radians [21,28]. 

Identified hadron acceptance and tracking efficiency were studied through 
Monte Carlo GEANT simulations [26,28,29]. At high p T (p T > 2.5 GeV/c) 
the efficiencies range from 73% to 87% and are nearly independent of px, but 
have a weak centrality dependence. Weak-decay feed-down (e.g. Kg — > it + it~) 
contributions to the pion spectra were calculated using measured K s and 
A yields [6] and a GEANT simulation. The feed-down contributions to the 
pion spectra were found to be ~ 12% at px — 0.35 GeV/c and decreasing to 
~ 5% for p T > 1 GeV/c. The final pion spectra are presented after subtract- 
ing these contributions. The inclusive p and p yields are presented without 
hyperon feed-down corrections to reflect total baryon production. The correc- 
tions range from < 20% for p+p and d+Au data [26,28,29] rising to ~ 40% 
for central Au+Au up to intermediate px, and are estimated to be less than 
20% at high p T [21]. 

Systematic errors for the TPC measurements were particle type and px de- 
pendent. They include: uncertainties in efficiency (~ 8%); dE/dx position and 
width (10-20%); background from decay feed-down, ghost tracks and PID con- 
tamination at high p T (8-14%); momentum distortion due to charge build-up 
in the TPC volume (0-10%); the distortion of the measured spectra due to 
momentum resolution (0-5%). The systematic errors are added in quadrature. 
Systematic errors for the TOF data for 7r^, p and p spectra are similar at both 
energies and are about 8% [28,30]. The total systematic errors for 7r ± yields 
at both energies are estimated to be < 15%, and those for p and p are < 25% 
over the entire px range studied [16]. 
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Fig. 1. Midrapidity (\y\ < 0.5) transverse momentum spectra for n , p and p for 
various event centrality classes for Au+Au at ^/snn = 62.4 GeV. Also shown to 
study the energy dependence are the central 0-12% -k^, p and p spectra for Au+Au 
at y / SNN = 200 GeV. The insets show -ir~ /ir + at ^/snn = 62.4 GeV and p/p at 
y/s^ = 62 A (0-10%) and 200 GeV (0-12%). The errors shown are statistical, and 
the shaded bands reflect the systematic errors. 

3 Transverse momentum distribution 



Figure 1 shows 7r ± , p and p yields for Au+Au at 62.4 GeV for 0.2 < p^ < 12 GeV/ c 
and various collision centralities. The hadron spectra at highpx (pt > 6 GeV/c) 
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for ^/snn = 62.4 GeV are steeper than the corresponding spectra for ^snn = 
200 GeV; comparisons of central collision spectra at both energies are shown 
in Fig. 1. This steepness mostly reflects the difference in initial jet produc- 
tion at the two collision energies. For high-energy p+p and d+Au collisions, 
particle production at midrapidity is found to follow m T (= \Jp\ + mass 2 ) 
scaling [16,17]. Such scaling implies that initial parton distributions dominate 
the particle production process [31]. The possibility of mi-scaling in heavy-ion 
collisions has been discussed in Ref. [31]. However, such mx-scaling is not ob- 
served in the data for v /snn = 62.4 and 200 GeV Au+Au. This will be evident 
from p(p)/n ratios presented later. The absence of m T -scaling may reflect a 
modification of the initial distributions through both partonic and hadronic 
final state interactions at RHIC energies. 

At ^/snn = 62.4 GeV, ir~ /n + = 1.01±0.02 (stat), independent of p T within ex- 
perimental uncertainties (inset of Fig. 1) and collision centrality (not shown). 
Similar features were observed at 200 GeV [21]. The p/p ratios show a slight 
decrease with px (inset of Fig. 1) and are independent of centrality. The de- 
creasing trend is more pronounced at 62.4 GeV [21]. For p*? < 3 GeV/c, 
p/p = 0.44 ± 0.01 and 0.77 ± 0.02 at v /5i^ = 62.4 and 200 GeV, respectively. 
For p T > 6 GeV/c, p/p = 0.29 ± 0.02 and 0.70 ± 0.05 at v /5^ = 62.4 and 200 
GeV, respectively. 



4 Nuclear modification factor 



The nuclear modification factor is defined relative to peripheral collisions 
(.Rcp) or relative to nucleon-nucleon collisions (Raa) [2]: 

[d 2 N/p T dydp T /(N hin )Y entral 

Rcp{Pt) = 



[d 2 N/p T dydp T /{N Un )}P eri P heral ' 



where (Nun) is the average number of binary nucleon-nucleon collisions per 
event, and 

d 2 N AA /dydp T /(N hiQ ) 



Raa(pt) 



d 2 a pp /dydp T /a : 



inel 
PP 



The a™ 1 are taken to be 36 mb and 42 mb for ^s^ = 62.4 GeV and 200 
GeV, respectively [32]. The d 2 a pp / dydp^ at 200 GeV are from STAR mea- 
surements [16]; for 62.4 GeV we use a parametrization of ISR data [33] in 
which the n invariant yield for p+p at ^/snn = 62.4 GeV is parameterized as 
Ed 3 a pp ^ x /d 3 p = A(e a - p T +b -PT +p T /p )~ n , with A = 265.1 mb GeV- 2 c 3 ,a = 
-0.0129 GeV" 2 c 2 , b = 0.04975 GeV^c, p = 2.639 GeV/c, and n = 17.95. The 
uncertainty in yields associated with this parametrization is ~ 25%. 

Figure 2 (upper panels) shows the pr, centrality and v /snn dependence of Rcp 
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Fig. 2. Upper panels: Centrality and pt dependence of Rqp for tt + +tt~ and p+p for 
Au+Au at y^SNN = 62.4 GeV. For studying the energy dependence, the correspond- 
ing i?cp for central 0-12% Au+Au at ^/snn = 200 GeV are shown. Lower panels 
: R AA for vr++vr- at 62.4 GeV (0-10%) and 200 GeV (0-12%) compared to three 
model predictions [8,9,11] (see text for details). A 25% uncertainty is associated 
with d 2 a pp /dydp T at 62.4 GeV. The error bars are statistical; the shaded bands are 
the systematic errors. The systematic errors for the 20-40% centrality data are of 
similar order as those shown for the 0-10% data. The shaded band around Rqp = 1 
at px = 10 GeV/c in the top right panel reflects the uncertainty in (Abi n ) calculation 
for 0-10% collision centrality. 



for 7r + +7r~ and p+p for Au+Au. The bottom panels show the 7r + +7r~ -Raa for 
the 0-10% and 0-12% centralities at a/snn = 62.4 and 200 GeV, respectively. 
For a given energy there is a distinct difference in the pr dependence between 
the i?cp for 7r + +7r~ and the Rqy> for P+P at intermediate pr- The Rev for P+P 
has a steeper fall with pr compared to n + +it~ . At high pr the Rqp values 
are similar for baryons and mesons at both energies. The relevance of these 
measurements for understanding the energy loss of quarks, gluons and their 
interaction with the medium will be discussed together with the p/ir + and 
P/ti~ ratios in the next section. When compared as a function of centrality, a 
dependence is observed for Rqp for both 7r + +7r~ and p+p at a/snn = 62.4 GeV. 
It is found to be stronger for 7r + +7r~. A similar decrease in Rep values with 
increasing collision centrality was observed at ^/snn = 200 GeV [21]. The 
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Rcp values at ^/snn = 62.4 GeV are higher than at ^/snn = 200 GeV for 
Px < 7 GeV/c; beyond this px they approach each other; this feature may 
be due to the interplay of initial jet production and the gluon density. For a 
smaller initial gluon density at the lower energy, the Rqp values at the two 
beam energies may approach each other at high px due to a steeper initial jet 
spectrum at 62.4 GeV [18]. 

The charged pion i?AA (left bottom panel of Fig. 2) for 3.0 < px < 8.0 GeV/c 
at ^/snn = 62.4 GeV decreases with p T and approaches ~0.35 at p T = 8 
GeV/c. In contrast the Raa values at ^/snn = 200 GeV are fairly constant for 
px > 4.0 GeV/c (bottom right panel). The difference in the px dependence of 
Raa at the two beam energies is influenced by the energy dependence of the 
following: the shape of the initial jet spectrum, the parton energy loss, and the 
relative contributions of quark and gluon jets. The steeper fall in .Raa with p T 
at 62.4 GeV may be due to the steeper initial jet spectrum. The constant value 
of Raa at high px for 200 GeV indicates that the effect due to the shape of the 
jet spectrum seems to be compensated by the parton energy loss. In addition, 
as quarks are expected to lose less energy than gluons in the medium [18,19], 
a higher contribution of quark jets at 62.4 GeV compared to 200 GeV for the 
same p T (x T 2,4 /^x° ~ 3) may also cause a difference in the energy dependence 
of Raa versus px- The differences in the high-px dependence of -Raa at the 
two collision energies rules out xx-scaling for Au+Au [34,35] in contrast to 
the observations for p+p [16]. This is expected, as various additional non- 
perturbative and perturbative processes for particle production in heavy-ion 
collisions have distinct px and ^/snn dependencies. 

In Fig. 2 the charged pion Raa are compared to model predictions at both 
energies to study their dependence on the initial gluon density, the lifetime 
of dense matter and the mechanism of energy loss. The predictions shown 
do not agree with the data in the region 2 < p T < 4 GeV/c, indicating that 
non-perturbative processes may dominate hadron production in this px range. 
The dashed curves are from a set of calculations which are sensitive to the 
choice of initial gluon density [8,36]. Comparison at high px shows that the 
initial gluon densities (dN 9 /dy) are about 650-800 and 800-1175 from these 
calculations for Au+Au at ^/snn = 62.4 and 200 GeV, respectively. The lower 
dashed curves are for higher gluon density. In addition, theoretical studies 
also suggest that for a given initial density, the Raa{pt) values are sensitive 
to the lifetime (r) of dense matter formed in heavy-ion collisions [9]. The 
solid curves are predictions from Ref. [9] at ^/snn = 62.4 and 200 GeV with 
r = 10 fm/c (i.e. larger than the typical system size of ~ 6-7 fm). For 62.4 
GeV, also shown is a prediction with r = 4 fm/c (upper solid line). The 
comparison at high px shows that, for this model, the lifetime of the dense 
matter formed in Au+Au collisions at ^/snn — 62.4 and 200 GeV is comparable 
or larger than the system size. Further insight to the mechanism of energy loss 
is obtained by comparing the data to theoretical predictions (dotted curves) 
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Fig. 3. Raa(pt > 6 GeV/c) versus N pilTt for charged pions for Au+Au at 62.4 GeV 
and 200 GeV, and their ratio. The error bars are statistical and systematic errors 
added in quadrature. The solid curves are results of calculations with radiative en- 
ergy loss for two different initial gluon densities in central collisions at both energies 
and then following the dependence of ln(i?AA) °c [13]. The dotted curves are 

theoretical calculations based on a parton model (see text for details) [37]. 

of Raa from models that consider only collisional energy loss [11]. For ^/snn 
= 200 GeV, the model predictions of Raa at high p T are close to the measured 
values and similar to corresponding Raa values from models based on only a 
radiative mechanism for parton energy loss. However, collisional energy loss 
model overpredicts the experimental Raa values at a/snn = 62.4 GeV and 
shows a stronger dependence on beam energy compared to models based on 
the radiative process of parton energy loss. 

The centrality dependence of Raa at high £>t may provide information on the 
path length dependence of parton energy loss in heavy-ion collisions. Figure 3 
shows Raa (pt > 6 GeV/c) as a function of iVp art for 7r + +7r~ for Au+Au 
at 62.4 and 200 GeV. The Raa values decrease with N paTt at both energies. 
The data are compared to results of two types of model calculations. The 
solid curves are from a model that uses a radiative energy loss mechanism for 
partons propagating through the medium formed in heavy-ion collisions [13]. 
The model assumes the parton production cross section to be a power law 
type and a parton energy loss that depends on the initial gluon density, the 
path length traveled by the parton, and the transverse area of the region of the 
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collision. Such a model predicts the centrality dependence of high pr Raa at 

2/3 

a given beam energy to be of the form ln(i?AA) ~ N' t [13]. The calculations 
are done for a set of two different gluon densities at both energies. The data 
follow the predicted dependence at both energies down to low values of N paxt . 

The dotted curves in Fig. 3 are results from a pQCD based parton model 
in which parton interactions with the medium formed in heavy-ion collisions 
are reflected through the modification of its fragmentation function [37]. The 
partons in the medium lose energy by induced gluon radiation. In such mod- 
els, the parton energy loss depends upon the local gluon density and the total 
distance of parton propagation [37]. These predictions are in reasonable agree- 
ment with the data for most of the centrality classes studied. 

The difference between the two models becomes clearer when we compare the 
ratio of Raa (pt > 6 GeV/c) values at 62.4 and 200 GeV with the ratios from 
data. This is shown in the bottom panel of Fig. 3. For the data, the Raa versus 
-^Vpart at 200 GeV is first parametrized by a polynomial function and then the 
ratios of _Raa(62.4)/.Raa(200) are calculated. The pQCD based parton model 
overpredicts the measurements for the most central collisions. For the most 
peripheral collisions measured, the model calculations from Ref. [13] slightly 
overpredict the data. It will be interesting to compare the -Raa versus iVp art for 
models with collisional energy loss and see if they provide further constraint 
on mechanism of energy loss of partons in heavy-ion collisions. 



5 Baryon-to-meson and anti-baryon-to-baryon ratios 

Figure 4 shows the p/n + and p/ir~~ ratios versus pr for Au+Au 0-10% and 
0-12% (upper panels), and 40-80% (lower panels) centralities at v /snn = 62.4 
GeV and 200 GeV, together with theoretical predictions to be discussed. 

The fact that for central collisions the p/tt + and p/ir~ ratios are close to unity 
in the intermediate pr region at 200 GeV has been attributed to either quark 
coalescence [7,22] or novel baryon transport dynamics based on topological 
gluon field configurations [38]. The quark coalescence models predict a specific 
energy dependence for p/V + , being higher at ^/snn = 62.4 GeV than at 200 
GeV in the intermediate pr region; the energy dependence is reversed for 
p/ir~ [22]. On the other hand, the baryon junction model predicts a decrease 
in the ratio at intermediate pr with decreasing collision centrality at a given 
^/snn [38]. 

As Fig. 4 shows, at a given pr the p/ir + ratio for Au+Au at 62.4 GeV is 
larger than the value at 200 GeV in the intermediate pr range, whereas for 
p/ir~ the reverse occurs. This specific energy dependence of the baryon-to- 
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Fig. 4. p/tt + and p/ir~ ratios versus pr at ^/snn = 62.4 and 200 GeV for central 
(upper panels) and peripheral (lower panels) collisions. For clarity of presentation, 
the systematic errors (shaded bands) are shown for only one of the beam energy for 
a given ratio. They are of similar magnitude at the other beam energy. The curves 
are model results [22,38,39] and are discussed in the text. 

meson ratio as a function of pt is consistent with the general expectation 
from quark coalescence models [22] . Our results also show that the baryon-to- 
meson ratios, here p/vr + , for the region 1.5 < pt < 6 GeV/c are higher than 
in p+p and d+Au [16]. This enhancement increases with centrality for both 
beam energies. 

The ratios for the 0-10% and 0-12% centrality data (upper panels of Fig. 4) 
are compared to predictions from models based on quark coalescence and a 
jet fragmentation mechanism for particle production at 62.4 GeV [22] and 200 
GeV [39], and baryon junction and jet fragmentation at 200 GeV [38]. For 
the intermediate px region there is a lack of quantitative agreement between 
model results and data. The recombination models predict a shift in the peak 
position of the ratios to higher p^ at the 62.4 GeV, which is not observed. The 
p/n~ ratios for the two energies do not cross-over as predicted by the models. 
The baryon junction model predictions are not in quantitative agreement with 
our 200 GeV data. 

At higher px the p/ir + and p/tt~ ratios are nearly independent of centrality at 
both 62.4 and 200 GeV. This observation, taken together with a constant Rep 
beyond p^ > 6 GeV/c, may reflect the dominance of particle production from 
the fragmentation mechanism. Also, at high p^ {pt > 6 GeV/c) we observed 
a similar p/ir~ ratio in central Au+Au and d+Au at 200 GeV [21,16] and a 
similar Rep for p+p and n + + 7r~ for Au+Au (see previous section). These 
observations appear to be inconsistent with the naive expectations from the 
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Fig. 5. The p/p ratios versus p T at y/s^ = 62 A (0-10%) and 200 GeV (0-12%). 
The errors shown are statistical, and the shaded bands reflect the systematic errors. 
Model predictions are shown as solid and dashed curves for 200 GeV [38] and 62.4 
GeV [22] central Au+Au, respectively. 

color charge dependence of the parton energy loss [18,19]. The difference in 
quark and gluon energy loss would have led to a lower p/tc~ ratio for Au+Au 
at high pt than that for <i+Au collisions and a lower Rep for p+p compared to 
7r + + 7r~. Recent theoretical calculations suggest that a much larger net quark 
to gluon jet conversion rate in the QGP medium is needed than given by the 
lowest order QCD calculations to explain the high p-£ particle ratios [20] . 

As one can see in Fig. 4, the jet fragmentation prediction is reasonable at high 
Pt for the p/ir + ratios at 62.4 GeV. However these calculations predict a much 
lower value for the ratio at 200 GeV. The failure of these model calculations at 
high pt is further noticeable when we compare them to the measured p/p ratios 
at both energies. Figure 5 shows the p/p ratios versus pt at 62.4 and 200 GeV. 
The data are compared to a model result in which baryons and anti-baryons 
are produced through baryon junctions and jet fragmentation at 200 GeV [38] 
and through coalescence and jet fragmentation processes at 62.4 GeV [22]. 
Both the models overpredict the data at lower (pt < 5 GeV/c). For p^ > 6 
GeV/ c, where fragmentation is the dominant mechanism of particle production 
in the models, they underpredict the measured p/p ratios at the two beam 
energies. The model calculations do not use the recent fragmentation functions 
for p+p as supported by the RHIC data from p+p and d+Au collisions at 200 



To further investigate the energy dependence of baryon-to-meson ratios, we 
present the ratio of p/rr~ between 62.4 GeV and 200 GeV and the ratio of 
p/tt + between 62.4 GeV and 200 GeV. Figure 6 shows that this double ratio 
of p/7i~ is independent of with a value around 0.6 for p T > 2 GeV/c, while 
the double ratio of p/ir + is around 1.2 for pt — 2 — 5 GeV/ c and increases with 
Pt, possibly due to different valence quark contributions at the two energies. 



GeV [16]. 



15 



2.5 



p/ji + 0-10% 
p/;t + 20-40% 
pin* 40-80% 



2t 




° p/jr- 0-10% 
Or p/Jt" 20-40% 
I A p/r 40-80% 



1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 



1 2 3 4 5 6 7 



8 



Transverse Momentum p T (GeV/c) 



Fig. 6. The ratios p/i: + at ^snn = 62.4 GeV and 200 GeV and p/vr~ at ^saw 
= 62.4 GeV and 200 GeV as functions of for various collision centrality classes. 
The error bars are statistical and shaded bands are the systematic errors. 



Baryon and meson production at high px and the relative contributions from 
quark and gluon jets have been discussed in Refs. [21,16]. The new observations 
presented in this paper necessitate further understanding of the role of gluons, 
quarks and their energy loss mechanisms. Gluon jets tend to produce more 
baryons than quark jets, whereas quark jets contributes substantially to pion 
production. This feature is supported by a much lower p/tt~ ratio at high 
Pt compared to the p/ir + ratio for low-energy p+p and p+A collisions [16]. 
The double ratios in Fig. 6 are consistent with this picture as the p/n~ ratio 
is lower at 62 GeV than at 200 GeV. Due to their larger coupling gluons 
should lose more energy in the dense medium formed in heavy-ion collisions 
than quarks [18,19]. This would lead to a lower p/n~ ratio for central Au+Au 
relative to peripheral Au+Au at both beam energies. This is not observed for 
the data reported here. A larger value of the p/7r + ratio at 62 GeV than at 
200 GeV is observed. This may be due to greater valence quark contribution 
at the lower beam energy. However, the p/ir + double ratio shows no centrality 
dependence. This is not expected if valence quarks contribute significantly 
more at lower energy and lose energy in the dense medium formed for central 
Au+Au. 

At intermediate p^, the features of the double ratios are not expected from 
the coalescence model; as seen in Fig. 4 the quark coalescence models will lead 
to more prominent baryon enhancement at 62 GeV than at 200 GeV. It is, 
however, surprising that the scaling is independent of centrality and extends 
to high when baryons are more enhanced at intermediate pp- 



16 



6 Summary 



We have presented a study of the energy dependence of 7r , p and p production 
for Au+Au at ^/s^ = 62.4 and 200 GeV. The pt spectra are measured 
around midrapidity (|y| < 0.5) over the range 0.2 < pp < 12 GeV/c. These 
measurements provide new experimental data for investigating the production 
of quarks, gluons and their interactions with the medium formed in heavy-ion 
collisions and the interplay between coalescence of thermal partons and jet 
fragmentation. 

The pt dependence of Rqp for charged pions and for protons and anti-protons 
is different at both energies. However, at higher p T the values of R C p for 
baryons and mesons at both energies are similar. The comparison of Raa 
versus pp to model predictions provides important information on quantities 
like initial gluon density and lifetime of dense matter. 

The p/ir + ratios for Au+Au at ^snn = 62.4 GeV are higher than the corre- 
sponding values at ^/snn = 200 GeV in the intermediate pr range, but the 
P/ti~ ratios are smaller. There is serious quantitative disagreement between 
data and the available theoretical models. We observe a scaling of the p/n~ 
ratios between corresponding centralities for the two beam energies at pt > 2 
GeV/c despite the strong centrality and p T dependence of these ratios. 
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